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Complete tensor-scalar and hydrodynamic equations are presented and integrated, for a self- 
gravitating perfect fluid. The initial conditions describe unstable-equilibrium neutron star configu- 
ration, with a polytropic equation of state. They are necessary in order to follow the gravitational 
collapse (including full hydrodynamics) of this star toward a black hole and to study the resulting 
scalar gravitational wave. The amplitude of this wave, as well as the radiated energy dramatically 
increase above some critical value of the parameter of the coupling function, due to the spontaneous 
scalarization, an effect not present in Brans-Dicke theory. In most cases, the pressure of the collaps- 
ing fluid does not have a significant impact on the resulting signal. These kind of sources are not 
likely to be observed by future laser interferometric detectors (such as VIRGO or LIGO) of gravi- 
tational waves, if they are located at more than a few 100 kpc. However, spontaneous scalarization 
could be constrained if such a gravitational collapse is detected by its quadrupolar gravitational 
signal, since this latter is quite lower than the monopolar one. 
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I. INTRODUCTION 



In order to test general relativity, one has to compare it to other, alternate theories of gravitation. Tensor-scalar 
theories, in which gravity is described by a spin-2 field combined with one or several spin-0 fields, are not only alternate 
theories, but "generalize general relativity" (see ||l| ); meaning that general relativity is obtained in them by setting 
all scalar fields to zero. Several such theories have been developed, from Fierz Jordan Brans and Dicke Q to 
Bergman [||, Nordtvedt Wagoner Q and, more recently, Damour and Esposito-Farese^jl] . In all these theories, 
the spin-0 and spin-2 fields {(p and 5^jy) are coupled to matter via an effective metric tensor = a'^{ip)g^^. The 
Jordan-Fierz-Brans-Dicke theory has only one free parameter whereas for Bergman, Nordtvedt and Wagoner the 
parameter is a function w(iy9). Damour and Esposito-Farese considered an arbitrary number of scalar fields, coupled 
one to the other. All these theories are motivated by, mainly, two theoretical reasons: 

1. they represent the low-energy limit of superstring theories ( |8| and [^) 

2. they give rise to new "extended" inflationary models pO[ | 

Since in Brans-Dicke theory \og[a{Lp)) is a linear function of 1^9, solar system experiments (weak field) are sufficient 
to constrain the theory, even in strong fields. Nevertheless, a more general theory, in which log(a((/j)) is a parabolic 
function (depending on two parameters), shows non-perturbative effects in strong field [ pl| , described as "spontaneous 
scalarization" in |12| . Thus, when describing neutron stars, general relativity and tensor-scalar theory can give 
significant differences for their masses, radii and gravitational field, whereas the difference can still be negligible in 
our solar system data. As a consequence, weak-field experiments cannot give much information on strong-field regime, 
and one needs to test this strong-field regime by other means: 

_ by looking for the orbital decay of binary-pulsar systems. This has been done by Damour and Esposito-Farese 
( ||l| and ||l^), who constrained the parameter space of the coupling function, 

_ by looking for monopolar gravitational radiation from collapsing compact sources which could be detected by 
the laser interferometric gravitational wave observatories (such as VIRGO [|l^ and LIGO |p^ ). 

This latter method requires that the signal be known and that the observed (or unobserved) amplitude be related 
to the coupling function parameters. Such kind of computations have already been performed by various groups, but 
they all considered only an Oppenheimer-Snyder collapse (i.e. "dust" matter, with no pressure), either in Brans-Dicke 
theory ( |l^ and |l^) or by doing some Taylor expansion of the coupling function ||l^. In this latter, the parameter 
space of the coupling function was restricted to the part where non-perturbative strong-fields effects do not happen. 
The aim of this paper is to present the results of computations of a spherically symmetric collapse, of a neutron star 
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toward a black hole Q with one scalar field and an arbitrary coupling function. All the hydrodynamics and field 
equations are treated with no approximation in order to get the monopolar gravitational wave form and amplitude. 
Moreover, including the equation of state allows t o star t the collapse with quite a realistic neutron star configuration 
and thus, spurious waves signals are avoided (see |III C| ). 

The paper in organized as follows. Section || describes the evolution equations for the star; section ^ gives the 
numeri cal re sults: initial- value models ( III A|) , collaps e and resulting wave signal ( III E ), comparison with previous 
works (IIIC) and exploration of the parameter space (HIE); finally section IV gives some concluding remarks. 



II. FIELD AND HYDRODYNAMIC EQUATIONS 
A. General equations 

As it has been stated before, the most general theory containing a spin-2 field and one (massless) spin-0 field 
contains one arbitrary coupling function a{ip). The action is given by 



S = (167rG, 



(2.1) 



where all quantities with asterisks are related to the "Einstein metric" g*^: G, is the bare gravitational coupling 
constant, i?* — g*'^ R*^i, the curvature scalar for this metric and 5, = det(f;^i,). The term Sm denotes the action of the 
matter, represented by the fields which is coupled to the "Jordan-Fierz" metric g^i, — a'^{ip)g'^^; all quantities with 
a tilde are related to this metric. That means that all non-gravitational experiments measure this metric, although 
the field equations of the theory are better formulated in the Einstein metric. The indices of Einstein frame quantities 
are moved through Einstein metric, whereas those of Jordan-Fierz quantities are moved through Jordan-Fierz one. 
By varying 5, one obtains 



R* 



2df,ipd„ip - g*g'i''dp(pda(p + 



SttG 



(2.2) 
(2.3) 



where 



a{ip) 



dS„ 



9* "yfii' 
d\na{ip) 



dip 



(2.4) 
(2.5) 



and Dg, = V*V* is the Laplace-Beltrami operator of g*^, V* denoting the Levi-Civita connection of 5*^. One 
can see that a{ip) is the basic, field-dependent coupling function between matter and scalar field. General relativity 
is obtained for a{(p) 0. 

The physical stress-energy tensor T'"' = 2{—g)^^/^ 6Sm/5g^i, is related to the Einstein-frame one by 

Tit = a\^)fl^ (2.6) 

The equations of motion are given by the stress-energy balance equation, written in the Jordan-Fierz frame 

V,f; = (2.7) 
and in the Einstein frame VlTi^"" = a(^)T*V5f(^ (2.8) 

Finally, let us call ipQ the cosmological value of the scalar field, which enters the theory as the boundary condition on 
the scalar field at spatial infinity. 



^numerically, the black h ole is never obtained, but the monopolar gravitational waves, far from the source, behave as if the 
black hole had formed (see [II E| ) 
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B. Coordinates and variables 



The present calculations have essentially been done by generalizing a previous work by Gourgoulhon |18| to tensor- 
scalar theory. Therefore, only a very brief presentation of coordinate and variable choice will be given here. The 
Einstein-scalar equations have been decomposed in the 3-1-1 formalism 1 19 onto a family of spacelike hypersurfaces 
labeled by the real index t called the coordinate time. The polar time slicing has been chosen in order to have 
good singularity avoidance (see e.g. |20| for discussion). On each hypersurface the radial gauge has been chosen 
with spherical-like coordinates (r, 9, </)), since the considered problem is spherically symmetric. All these assumptions 
(spherical symmetry, radial gauge and polar slicing, called RGPS) imply that the metric y*^, = a^"^ {ip)g (which 
verifies Einstein-like equations (2.2), with an extra-term) takes the diagonal form: 

ds^ = ~N'^{r,t)dt^ + A^{r,t)dr'^ +r'^{d9'^ +sin^ 9d(j)'^) (2.9) 

where N{r, t) is called the lapse function. The metric g*^, will often be described by the three functions i^(r, t), m{r, t) 
and C(r, t) defined by 

A^(r,t) = exp(zy(r,t)) (2.10) 

, / 2TO(r,t)\"^/^ 
A(r,t)= 1 ^ (2.11) 



and C(r,t) = In (2.12) 

All coordinates are expressed in the Einstein-frame, and asterisks are omitted. However, "physical" quantities will 
often be written in the Fierz metric and noted with a tilde. 

In this work, neutron stars are modeled as self-gravitating perfect fluids. They can be considered to be made of 
degenerate matter at equilibrium, the equation of state being temperature independent (cold matter). This docs not 
hold only soon after their formation. The stress-energy tensor writes 

T^u = (e + p)Uf_iUi, + pg^iy (2-13) 

where is the 4- velocity of the fluid, e is the total energy density (including rest mass) in the fluid frame and p is the 
pressure. The relation to its Einstein-frame counterpart is TIf — a'^{Lp)Tlf. The description of the fluid is completed 
by an equation of state 

e = i{fiB) (2.14) 

with fiB being the baryonic density in the fluid frame. One then deduces the pressure as a function of n^. Let 
r ~ NvP be the Lorentz factor connecting the fluid frame and Et hypersurface frame, calling 

E ^ one gets E ^r^{i + p) ~p (2.15) 

The fluid baryonic number is represented by the coordinate baryonic density 

~ number of baryons in 6V 

D = — = ATuB (2.16) 

oV 

where SV = r^ sin 9drd9d(l) is the element of the coordinate 3- volume on a given Ei, defined as the set of points whose 
coordinates are between r and r + dr, 9 and 9 + d9, (j) and (/) + d(j). The fluid motion is described by the following 
variables 

dr 

V = — — — (coordinate velocity) (2-17) 

dt v}' 

proper distance travelled on Et A , ^ 

U^—T \ ■■ ^T-^^irr^ 2.18 

elapsed proper time on Ej N 



one then has F = (1 — C/^) and deduces the components of T^n given by (2.12). The fluid log-enthalpy is also 
introduced, defined as: 

H = \n( /^^^ ) (2.19) 
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and, finally, three "scalar-field" variables: 



N~dt 



^ = 77 



(2.20) 

(2.21) 
(2.22) 



C. Tensor-scalar field equations 



Spherical symmetry helps t o obtai n gravitational field equations; we followed the procedure described by Gour- 
goulhon [Q, projecting eqs. (2.2-^) on the 3-surfaces St and along their normal. Hereafter, we use the following 
notation: 



SttG, 



The tensor Einstein-like equations (2_^) turn then into: 
_ one Hamiltonian constraint equation 

dm 



(2.23) 



three momentum constraint equations which reduce to only one non-vanishing 



dm 
'dt 



2^ij7^-q^a^{ip){E + p)V 



(2.24) 



six Einstein dynamical equations which here reduce to two non-vanishing; one is degenerate, only giving a 
condition on the lapse function 



dv QttA'^ 
Ih ~ 2 



Aur' 



+ a-^{ip)r{p + U'^{E+p))+rE 



(2.25) 



The other one will not be used in this work. 
Writing the scalar-field wave equation ( |2.3| ) with our variables gives 



(2.26) 



One more equation concerning the scalar field will be used; although it is redundant with (2.23)-( 2.26| ), from which 
it is deduced , it will be useful for numerical integration: 



2~dt 



N 
~A 



or 



^i^V^ - \ + i"lM^V{v)N{E -3p-{E + p)U^) (2.27) 



D. Matter evolution equations 



In order to get the evolution of the variables E and [/, let us consider the momentum-energy conservation equation 
(2.8), we get 
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dU_ 
'dt 



dE 

dt dr 
dU 1 



+ — — {r^{E + p)V) =-{E + p) {a{ip)N [(3 + U^)i; + AUtj] + rAN [(1 + U^)^tj + UE] } 



V 



dr 



E + p 



U 



dp N dp 
dt ^ Id^ 



AN 



A 



2^ 



(2.28) 
(2.29) 



Equation (2.28) expresses the total energy conservation (matter plus gravitational and scalar energy), eq. (2.29) being 



the tensor-scalar analogous of the Euler equation. One notes, in the latter, the -^^^ + ^If) term which may 

cause some trouble when numerically calculati ng it, since it is the quotient of two quantities vanishing at the surface 
of the star. Thus, if one uses the log-enthalpy (|2.19|) , this term may be replaced by j4- + T^)' "^tiich is well 

defined near the surface. 

Expressing the baryonic number conservation 



VnUBU^ = 



(2.30) 



one obtains 



dD 
'dt 



+ aM^^ (r^DV^ + a{ip)DN{3ijj + AUr]) = 



(2.31) 



Because log-enthalpy is used in order to avoid numerical singularities at the surface, there has to be an evolution 
equation of that quantity. Since H = H{nB), and ub ~ one may write 



dH 
'dT 



dH dH fdfi 



dr 



dfiB V 
dH V 



dfiB AT 



(1 



V- 



dhi 



dr 



, dD dH 
dr driB 



1 fdA ^dA 

— h V — 

A\dt dr 



--If 



dr 



(2.32) 



with terms in the right-hand side being replaced using eqs. (2.23), (2.24) and (2.29) by source terms involving 
E^p and S. The results of numerical integration of all these equations (2.22)-(2.32) will be presented in next section. 



III. NUMERICAL RESULTS 



The numerical procedure, the code and its tests are described in Appendix ^ In this section, only results are 
presented and discussed. An important choice is that of the coupling function &((/?). Following p^ , we chose a 
function depending on two parameters, for all our study: 



a{if>) 



_ gao(¥'-¥'o) + %-(v-Vo)^ 



(3.1) 



Figure 9 of jp^ gives constraints on the (aoj/3o) space of parameters, imposed by binary-pulsar measurements. 
Section [II D investigates this space of parameters for scalar gravitational waves. Note that Brans-Dicke theory is 
obtained for (3q — 0. 



A. Static configurations 



Physical scenarios to form a black hole involve either an accreting neutron star or a post-supernova remnant (when 
a part of the ejected envelop falls back onto the new-born neutron star). In both cases, the mass of the neutron 
star must reach its maximal value above which the star becomes unstable. It is then interesting to get unstable 
equilibrium configurations of neutron stars, endowed with a scalar field, close to the maximal mass. They are used as 
initial configurations for the collapse. Thus, setting all d/dt terms to zero, as well as V and U in eqs. ( 2.23| )-(2.28), one 
gets the scalar equivalent of the Tolman-Oppenheimer-Volkoff system. The system obtained is the same as equations 
(7) in 1 11 1, since the same gauge is used. Considering a polytropic equation of state 
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e{nB) = fiBriiB + K — 

7-1 V «o 



TIB ' 



p = KfiQiTiB — (3.2) 

with ?7iB = 1.66 X 10^^^ kg and uq = 0.1 fm^'^, one can integrate the scalar TOV system, starting at the center with 
a given value for hB{r — 0), up to the surface at which hB{r — Rstar) — 0. In this study, two types of polytrope will 
be used: 

1. 7 = 2.34 and K = 0.0195, which has aheady been used by ||ll| to fit EOSII in ||l|l; called EOSl in this paper. 

2. 7 = 2 and K = 0.1, as used in ||2^; called E0S2 in this paper. 

Then for each static configuration, the total ADM mass of which will be called gravitational mass and the total 
scalar charge w such that, for r ^ oo, (p{r) = ipo + G,,uj/r + 0{l/r'^), can be determined through eqs.(8) of ||ll|. These 
two quantities are useful to match the obtained interior solution to the exterior one (spherically symmetric solution in 
vacuum), which is known analytically in an other gauge (described in p^]) and thus one can obtain a static solution 
everywhere. The resulting fields are shown in Fig. for Pq — —6, with large value of the scalar field inside the star, 
even for a very small asymptotic field cpo (spontaneous scalarization) . These solutions arc then used as initial values 
for the dynamical evolution. They have been computed, with increasing central densities ^^(r = 0), in order to get 
an "unstable" configuration (for which gravitational mass is a decreasing function of the density). This property is 
not evident in tensor-scalar theory, but the dynamical code being sensitive enough to trigger the instability only by 
round-off errors, it has be en c hecked numerically. The hydrostatic equilibrium is obtained, thanks to pseudo-spectral 
techniques (see Appendix A 1 and [p4[), up to very high accuracy (10^"'^° relative error on the hydrostatic equilibrium). 



which enables the dynamical code to be sensitive to instability (see ]lq]). 



B. Scalar gravitational waves 



Hereafter, four collapse calculations will be presented, called A, B, C and D. The parameters of the static config- 
urations, which were used as initial conditions for the collapses, are described in table |. Note that collapses A and 
B use a 7 = 2.34 polytrope, whereas C and D use a 7 = 2 one. First, only the case A will be considered. As far 
as the hydrodynamic part is concerned, the collapse is very similar to that in general relativity, described in fl^ . It 
can be seen from Fig. ^, that A(t, Rstar (t)) 00 due to the pathological behavior of the radial gauge when Rstar is 
approaching the Schwarzschild radius. An apparent horizon is expected to develop, but minimal 2-surfaces cannot be 
described by the radial gauge. On Fig. || are plotted several quantities during the collapse (until N{r — 0) becomes 
too small). Thus, although RGPS coordinates are not well adapted for the description of a black hole, they were used 
to describe the collapse toward it, as in Moreover, from Figs. ||and|^ one sees that the star has almost entered its 

Schwarzschild radius is — = 1.001 at the end of the collapse); so that no significant later evolution could be 

achieved inside the star. Actually, one notices that the lapse goes to zero within the Schwarzschild radius of the star. 
Since all evolution equations are written d/dt = N x (source term) and the coordinate velocity V — (N/A) U 0, all 
hydrodynamic and scalar-field quantities are "frozen" inside the star. Therefore, their evolution can be numerically 
stopped, in order to avoid the singularity of A{r = Rstar)- However, all field quantities continue to evolve outside the 
star, as long as one wants in terms of coordinate time (which is the time of an observer at spatial infinity). 

The results of this evolution are shown in Fig. |^, for collapse A. The fate of the scalar field is particularly interesting: 
for r 3> Rstar the field relaxes toward the asymptotic constant value set by cosmological evolution; the scalar energy 
of the star is radiated away as scalar gravitational wave. The scalar field (</?(r = Rout)) is considered to be sufficiently 
far away from the star (i.e. in the wave zone) to give the monopolar gravitational wave signal. Using the "frozen 
star" to evolve the fields outside the star, makes the integration time long enough to get all the information from the 
collapse at r 3> Rstar (where all gauges become equivalent). The main difference from general relativity is the scalar 
monopolar radiation, which carries away energy and can interact with a detector. Looking far from the source (at a 
distance r ^ Rstar), one can write the metric (see M) 



9ti.y{r-,t) = a^((po) 



ft,. + - (/im- + 2aof^/p.) + 0{r- 



(3.3) 



where /^^ is the flat metric, h^^{t — ^) and F{t — ^) are respectively the quadrupolar and monopolar components of 
the wave. Since this work is done in spherical symmetry, only the monopolar mode shall be considered. The function 
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F is the same as that of eq. (A3), and is plotted for collapse B in Fig. |^. If one wants to compare the amplitude to 
those of general relativistic gravitational waves (hfj^), at a distance d, then the right quantity is (see jl^ or p5|): 

(3.4) 



h{t) = -a\ipo)aoF{t) 



The Fourier spectrum of the signal is 



Hf) 



(3.5) 



and the power spectrum {f\h\), which is plotted in Fig. |^ for the collapse B, is then useful to determine the characteristic 
frequency fc of the signal. Actually, the quantity to be compared with detectors' sensitivity is 



Kfc) 



(3.6) 



which is expressed in Hz From Fig. |[ one sees that a "collapse B" could be detected by VIRGO, provided it 
occurs within 300 kpc. Finally, the radiated Einstein- frame Bondi energy writes 



Escal — 



G 



* Jo 



dF 
~dt 



dt 



(3.7) 



which is the total energy radiated in gravitational waves. 



C. Comparison with previous works 



Because all previous works have studied an Oppenheimer-Snyder collapse, for comparison some runs were done 
putting p and H to zero. 

The first work by Shibata et al. |jl^ considered the dust collapse in the Brans-Dicke theory (/3o = 0). The authors 
used two types of initial conditions, called {A) and {B). The case {A) starts the collapse with (p — ipo, whereas (B) 
starts it with a quasi-static solution for (p. Making the same dust collapses, the same waveforms and amplitudes were 
obtained (their <I> is related to i^s by $ = 2 log {A{(p)) = 2ao{(p — (fio))- Taking the equation of state into account, with 
equilibrium initial configurations (which are the most realistic possible) , gives the waveform of Fig. ^ (/3o = 0) . The 
form and amplitude are very close to the {B) type collapses of Shibata et al.(see their Fig. 3), the power spectra being 
close too (Fig. |], with lower amphtude, and Fig. 6 of |l5|]). 

Another work in Brans-Dicke theory, with dust matter, was done by Scheel et al. |]l6|, dealing more with the fate 
of the final object. However, comparing 2aodF/dt of this work, with /' of their paper gives the same result. Whereas 
the authors can conclude on the final state of the collapse, this work cannot describe it. All that can be said is that 
the scalar charge (or mass) is all radiated away, so that the scalar field relaxes toward its cosmological value, and that 
an apparent horizon is expected to appear (Cf. sec. HI B ). These results still hold in more general tensor-scalar theory. 
Finally, it can be pointed out that, contrary to their work, here the tensor mass in the Einstein frame is considered, 
not in the Jordan-Fierz one. 

The last study of spherically symmetric collapse in tensor-scalar theory was done by Harada et al. They used 
the Oppenheimer-Snyder metric (in general relativity) as a background spacetime and did some Taylor expansion of 
tensors and equations in terms of scalar field coupling function parameters. Therefore, they used unrealistic initial 
conditions {ip = ifo) and could not study the cases of spontaneous scalarization. Waveforms resulting from our 
calculations for different Pq, from —10 to 50 are shown in Figs, ^to ^ with unstable eq uilibri um initial configurations 
and full hydrodynamics, ag = 3.16x 10~^ as in [0 and the mass is the maximal one (Cf.III A), the ratio Rstar/M ~ 4. 
The results are different, although showing the same tendancy as in for many (3q, due to the fact that Harada 
et al. took unrealistic initial conditions. Thus in their simulations, when the collapse begins, the scalar field on the 
one hand evolves to reach its quasi-equilibrium configuration, on the other hand, it feels the effects of the collapse. 
Their signal is then a superposition of these two effects: raise of the scalar field up to its equilibrium value, with one 
or several oscillations depending on /3o, then the fall of this field, due to the collapse of the matter. The signals for 
/3o > 20 have more important oscillations than other ones; this may be explained as follows. Inside the star, there 
are unstable modes which develop for these Po and when the ratio R/M becomes small, as it has been shown by 
Harada in for static configurations (see Fig. 5 of p3|). This is due to the fact that, near the star's center, E ~ 3p < 
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allowing for spontaneous scalarization to develop^. Here, one can see an effect of the pressure on the signal. These 
non-perturbative modes develop on a time-scale r 



(see eq. (5.3) of [^) with rg' being the free-fall time of the star. Thus, for Po > 20 the modes develop during the 
collapse, when R/M becomes small enough. This rise is fast enough to be (at least partly) seen before the star enters 
its Schwarzschild radius (see Figs. || and p|). 

D. Exploring the parameter space 

Thanks to the fact that it solves complete equations, the code presented in this paper is able to explore a larger 



part of the parameter space and give more "physical" results. In sec. |IIC|/3o has been varied from -10 to 50 and 
was fixed. The effect of spontaneous scalarization, for /3o — —5 and lower, was observed changing the amplitude, but 
not the shape of the wave. However, the two regimes (depending on the value of /3o) were studied separately, when 
varying ap- This latter has been taken between and 3 x 10~^, either in spontaneous scalarization (/3o = ~6), or 
without it {Po = —4). The results are shown in Fig. ^for he and in Fig. ^for Escai- One notices that for /3o = —6 
he (X ckQ and Egcai constant, whereas both c>c a§ for /3o = —4. The scalar field amplitude does (almost) not depend 
on ao in spontaneous scalarization, whereas it is directly proportional to it otherwise, like in Brans-Dicke theory. 



The equation of state has also been changed. Collapse C and D have been performed with E0S2 (see sect. HI A ). 
Results are shown in Fig. |l^ for collapse D and are similar to those of EOSl (collapses B), the scalar gravitational 
wave signal having the same shape, shifted since fc is higher (1 kHz for D, versus 800 Hz for B). This difference as 
well as that in amplitude can easily be explained by the change of mass of unstable-equilibrium configuration. Thus, 
varying the parameters of the polytropic equation of state is not of much interest. Better would be the use of a 
realistic equation of state (as in ||l^) but, for this work, a polytrope gives already good results. 

Finally, a few more runs were performed, with EOSl, varying 0q from 50 to —10 and, for each f3o, the maximal 
value of ao allowed by solar-system experiments (see e.g. Fig. 9 of |1^) was taken. More precisely, the lowest value of 

1 2 

al < 10"^ and < — x 10"^ (3.8) 

IpoI 

was taken. Results are shown in Figs|l^ and ^for Escai, and it can be seen that effects of spontaneous scalarization 



appear for Pq < —4.4 and that, for f3o > 20, one can see effects of instabilities of the scalar field (see section HIC). 
As far as the gravitational wave signal is concerned, the quantity plotted in Figs. |l^ and ^ is he, at 10 Mpc from 
the source. The characteristic frequency range is 700 Hz < /c < 900 Hz except for > 20, where fc ~ 2000 Hz. The 
maximal signal he — 5 x 10~^^ at 10 Mpc which is, at least, one order of magnitude lower than the best expected 



sensitivity of interferometric detectors currently under construction (see e.g. |27 ). However, one may compare the 
signal amplitude and energy to those of similar collapses in general relativity(2 and 3D stellar core collapse, see 
p8| and [g9[) and see that both are quite higher. Thus, if spontaneous scalarization effects are likely to occur, 
their gravitational signal should be more easily detected than the quadrupolar one from a collapsing source (aborted 
supernova or neutron star reaching its maximal mass by accretion). 

IV. CONCLUSIONS 

This work has been d one with very few approximations (Ax N set to a constant at the outer edge of the grid, see 



eq. (A6) and Appendix A3, evolution "frozen" when the lapse becomes too small); all the tensor-scalar equations, 
including hydrodynamics, were solved with high accuracy by means of spectral methods. Although the gauge choice 
does not allow any study of the state of the final object, the scalar gravitational wave signals can be obtained and 
the results compare well with previous (simplified) works. From that comparison it can be said that, whereas taking 
pressure into account, in most cases, does not have significant effects on the signal (which was to be demonstrated), it 
is the only way to get reliable initial conditions and to trigger the collapse in a "natural" way. Doing so helps to get 



^see the simplified model of 
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clean waveforms. On the other hand, since the complete set of tensor-scalar equations was solved, it was possible to 
study the effects of coupling function parameters. Mainly, one sees that the outgoing monopolar gravitational wave is 
very dependent on the coupling function, especially the /3o parameter. This is interesting because the ao parameter 
can be constrained by solar-system experiments, since it represents only a linear deviation from general relativity, 
whereas f3o cannot be really probed by that mean. Even if the signal from extra-galactic sources is not strong enough 
to be detected, it is higher than the quadrupolar one and involves more energy. It means that, if a quadrupolar wave 
signal is detected by VIRGO or LIGO, with no monopolar component, then the constrain on tensor-scalar theory will 
be quite strong. A future project of study is the supernova collapse and bounce in this framework since, in that case, 
one has electromagnetic and neutrino signals which make possible the use of even negative results of detection. 
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APPENDIX A: NUMERICAL PROCEDURE 



Hereafter, some details of the numerical techniques (spectral methods) are described. More complete explanations 
can be found in ||, |o| and 



1. Chebyshev decomposition 



The numerical problem is to solve a set of partial differential equations. For this purpose, each field (or function) 
/to(r) is represented, at a given time to, as a truncated series of Chebyshev polynomials (or as a column vector of the 

Usual number of coefficients (or points) is between 17 and 65. Then all spatial operators 



coefficients of this series), 
like 



^-^1^' / 

or I r 



reduce to matrix multiplications of the set of /'s coefficients. Constraint equations ( 2.23 ) and ( 2.25 ) are thus easily 
integrated. Other operations, like multiplication of two functions, must still be done in the physical space (at grid's 
points). 

Evolution equations arc written in the form 



dl 

dt 



fi. 7+1 ft. J 



/2 



(Al) 



with S being a spatial operator on / and tj being the J-th instant of integration. There is need to evaluate S at 
time ij+i/2, which can be done either explicitly (extrapolated from known quantities at tj and tj_i) or implicitly 
(interpolated from unknown quantity at ij+i). Explicit integration suffers from the severe Courant-Friedrich-Levy 
constraint on the time step. However, in the case of an advection equation (like (2.31)), it can be used almost with 
an arbitrary time step provided that the advection velocity is zero at the edges of the grid; hence the comoving grid 
with the fluid (see Appendix |A 3|) . The implicit way requires to write 



1 - 



J+i 



^ -/tj+i — ftj + 



J+i 



tj) 



Sift,) 



(A2) 



where the function / is represented by its coefficients and with (1 — ^''' S) being a matrix, which is inverted to 

get the solution ftj^i- In the case where one has to impose boundary conditions, the right-hand side of eq. (A2) is 
replaced by a vector of coefficients containing zeros except on the last or last but one column. One thus gets a "free" 
solution which can be combined with the first one in order to satisfy boundary conditions. 



-tj) 
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2. Wave equation 



Unfortunately, the S operator present in eq. (|Al|) i s not always linear (i.e. not represented by a matrix in coefficient 
space). That is the case for the wave equation ( 2.26| ), with a spatial operator of the form ip — > e^^A(^. This equation 
is then decomposed as follows: 



with 



(e^^ - Xi*\^ + V + Xl*^] Aip + e'^^^^ + + f Matter source terms') 

v / ar or at at \ I 



and A2 + Aj^ r + Aq being an approximation of e^^^''*', allowing to write the most important part of the spatial 
operator in a "linear" form. Finally, one writes d"^ ip j dt^\^_^ ^ with second-order approximation and Apt, — (^'/'tj+i + 
Apt ,_-^)/2 to make the integration scheme implicit. The boundary condition imposed on the outer edge of the grid 
(far away from the star) is that of an outgoing wave, meaning that the wave can be written as 



p{t,r)^^o + -Fit--) 
r c 

which is an exact condition for 1-D waves. Differentiating, one gets 



1 dp dp 
c dt dr 



(A3) 



(A4) 



outer edge 



Actually, this is not the right boundary condition for a wave equation in curved space-time, however, since the 
boundary condition is imposed far away from the star (i.e. on a nearly flat space-time), a good approximation is 
obtained by taking as a boundary condition 



-C 



dp 
'dt 



dp 
dr 







(A5) 



outer edge 



3. Dynamical evolution 



The integration procedure is quite similar to that of Gourgoulhon in [Q. All matter and field quantities are 
supposed to be known at some initial instant to, and one wants then to get t hem at to + 5t. First, one can compute 



the scalar field variables p and S at that time thanks to wave equation ( 2.26 ) and eq. ( 2.27 ). Similarly, one gets the 
fluid quantities E, U (and thus F), H and £>, with their evolu tion e quati ons (b.28D , ( |2.29| ), ( ^3^ and ( |2.3l| ). Then, 
one can deduce the metric coefficient A(r, tg + 5t) through eq. (2.23) and (2.11) (the integration constant is ob taine d 
by the condition Vt, A{r = 0,t) — 1); and determine the fluid proper baryonic density n^, by inverting relation (2.16). 
The equation of state t hen g ives the pressure pifiB)- 



Finally, one uses eq. (2.25) to obtain v{r,to -\- St), which is determined up to an additive constant. Since there is 
no Birkhoff Theorem in scalar-tensor gravity, this constant cannot be determined by matching the interior solution 
to the exterior (static) one, used before. Even the exterior space-time is dynamic. Using a large grid (which is going 
far away from the star, typically '--^ 30 x Rstar) enables to be in weak-field regime at the outer edge of the grid and 
so, to write with a good approximation 



yt,A{Rout,t) X N{Rout,t) = Kan 



(A6) 



Kan being a constant determined for the static configuration {Kan = 1 for General Relativity). One thus gets the 
integration constant for i^(r, to + St). Observing A{Rout) dur ing the collapse, one sees that AA{Rout) < 10~^. This is 
quite lower than the overall committed error, see Appendix A 4. Once N{r,t) is obtained, the velocity V is deduced 
and all quantities are known at the instant to + St. 

Time integration is performed by means of a second-order semi-implicit scheme (see previous sections), boundary 
conditions being imposed on the H, U system 



]_d_ 

J.2 Qj, 



r^V 







(A7) 



^Rsta 
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and on the wave equation (A5). Contrary to iQ, the code is working with several domains. Each domain is partly 
comoving with the fluid, meaning that inside the star, both edges of a domain are comoving, the grid's velocity at 
each point being (linearly) interpolated. The edge of the outer domain is at rest [Vgridir — Rout) — 0), and the 
grid's velocity between star's surface and this outer edge is interpolated at each point the same way. Typically, two 
domains have been used for the star's interior and three for the exterior, with about 65 points in each domain, which 
is considerably lower than the number of points used in finite difference schemes. A full run, with this multi-domain 
spectral methods, took about 20 minutes of CPU time on an Onyx Silicon Graphics workstation (with MIPS R4400, 
200 MHz processor), for ~ 20000 time steps. 



4. Tests 



All equations have been checked using a Mathematica algebraic code. More precisely, the Schwartz relation for 
equations (2.23) and (2.24) 



drdt dtdr 







has been computed and one then gets 



dip 
'dt 



dE_ 
'dt 



right-hand side of eq. 



which is consistent with ( 2.26| ) and ( 2.28 ). Finally, setting ip — gives the equations of general relativity as described 
in 

For static configurations, the following test was performed: considering the same equation of state as in |11| (EOSl 
in this work), the same coupling function (a(</?) — exp{~3ip^)) and asymptotic scalar field value {ipo = 0.0043), 
we obtained the same dependence for the effective scalar coupling constant (the ratio between the scalar and the 



gravitational energies) on the star's baryonic mass (Fig. 2 of | 
gravitational mass of neutron stars, when taking into account t 



ll|). Wc also observed an increase of the maximal 
le scalar field, as it has been showed in Fig. 1 of [O. 



The subroutine solving the wave equation has been checked by taking analytical solutions of simpler wave equations 



W{t,rfAp 



(A8) 



where W{t, r) is a given function, and verifying that the discrepancy between numerical solution and the analytical 
one goes down as the square of the integration time step (second order scheme). For example, 



which is solution for W{t, r)^ = cosh(r)^/2(2 -I- 1)"^ \n{t + 2), was numerically obtained at 10~^ relative accuracy, with 
a time step which is dt' = 1/100 of grid's radius, after 1000 iterations. The boundary condition has been checked by 
looking for the remaining energy of the wave in the grid, after the wave was supposed to be "gone out" . With the 
same time-step as above, and for the usual wave equation {Dip = 0), after 500 integrations, there remained 4.9 x 10~* 
of the initial energy, this amount going down as {dt'^)'^, for the energy is the "square" of the amplitude. As it has 
been seen in previous sections, a numerical grid partly comoving with the fluid was used. For this purpose, the wave 
equation have been adapted to such a grid and tested. 

Finally, during dynamical evolution, a good test of overall accuracy was made by comparing the energy density e 
given by the baryonic density through the equation of state, to that deduced from E (which is an evolved quantity). 
The same is possible for H. On a whole run this error always stayed < few x 10~^. The conservation of the baryonic 
number was verified with a relative accuracy better than 10~^. 
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FIG. 1. Density (ns), metric potentials {A and N) and scalar field (ip) profiles for a neutron star of 2AMq, for EOSl 

(7 = 2.34 and K = 0.0195 polytrope) and with a coupling function 0(9?) = exp(— 3^3^). The asymptotic scalar field value is 
tpo = 10~®. Star's radius Rstar = 13.1 km. 
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FIG. 2. Profiles of various quantities at difTerent values of t between and 4.64 ms, for collapse A. The fluid velocity U{r,t), 
measured by the hypersurface observer, is expressed in units of c and its evolution is downward, the extremity of each curve 
giving the position of the star's surface at the corresponding instant. The evolution for A{r,t) (metric potential), H(r, (scalar 
"energy") and H{r,t) (log-enthalpy) is upward, for ip{r,t) (scalar field) and N{r,t) (lapse) is downward. 
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t [ ms ] t [ ms ] 

FIG. 3. Evolutions of various quantities during the collapse A, as a function of the coordinate-time t. Rout = 300 km is the 
radius of the outer edge of the grid. (p{r,t) is the scalar field and the radius is the coordinate value for which hs, the baryon 
density is zero. N(r = 0) is the lapse at star's center and V{r = R) is the star's surface velocity. 



15 




t [ ms ] 



FIG. 4. Waveform of the emitted signal during collapse B. The plotted quantity is the function F{t) (see sec. IIIB ), measured 
at r = 300 km and expressed i n m eters. ip{r, t) is the scalar field and ipo its asymptotic value. To get the gravitational wave 
amplitude h, one must use eq. (3.4). 



Fourier power spectrum (o<„ = 2.5 x 10 , ^ ci= — 5) 
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FIG. 5. Fourier power spectrum of the scalar gravitational wave emitted during collapse B. h is the Fourier transform of 
the signal, measured at 10 Mpc and /, the frequency. 
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t [ms] 

FIG. 6. Same as Fig. ^ but for different collapses, ao = 3.16 x 10^^ for all of them, but (3o varies from —4 (upper curve, 
— ) to l3o = (lower curve, — ■ ■ ■ — ). 




t [ms] 

FIG. 7. Same as Fig. ^ but (3o varies from —10 (upper curve, — ■ — ) to /3o = —6 (lower curve 
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FIG. 9. Same as Fig. | but /3o = 40 {—) and /3o = 50 ( ). 
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Scalar wave amplitude at lOMpc 




«0 

FIG. 10. Amplitude of the scalar gravitational wave he (see eq. (3.6)), at 10 Mpc, as a function of the coupling coefficient 
ao, for two different values of /3o. 
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FIG. 11. Radiated scalar gravitational energy E^cai 
coefficient ao, for two difi'erent values of /3o- 



(see eq. (3.7)) emitted during a collapse, as a function of the coupling 
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power spectrum («q = 2.5 x 1 O , 13 — 5) 
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FIG. 12. Same as Fig. ^ but for collapse D. 




FIG. 13. Radiated scalar gravitational energy Escai (see eq. (3.7)) emitted during a collapse , as a function of the coupling 
coefficient /3o, between —6 and 0. For each /3o, Q-o has the maximal value, imposed by solar-system experiments (see eq. (3. J)). 
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Scalar wove amplitude at 1 OMpc 
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FIG. 15. Amplitude of the scalar gravitational wave he (see eq. (3.6)), at 10 Mpc, as a function of the coupling coefficient 
, between —6 and 0. For each /3o, ao has the maximal value, imposed by solar-system experiments (see eq. (3.8)). 
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Scalar wave amplitude at lOMpc 



in 




FIG. 16. Same as Fig. ^ but /3o varies from to 50. 
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TA BLE I . Initial condition parameters of the collapses presented in this paper. The equations of state (EOS) are described 
in sec. |lll A , tpo is the asymptotic scalar field value (given by cosmological evolution), ao and /3o are the coupling function param- 
eters (3.1), Rstar denotes star's radius, ns(r — 0)is the central baryon density (in units of nuclear density, 1 n„„c ~ 10*^ m-^), 
Mg is the p*,^-frame ADM mass, Mb the baryonic one and lo the scalar charge. 
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